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Abstract 

We obtained center conditions for a O-symmetric system of degree 5 
for which the origin is a uniformly isochronous singular point. 
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1. Let us consider a planar differential system 

x = y + xR n _i(x,y), 

V = ~x + yRn-i(x,y), 

where R n -\{x, y) is a polynomial in x, y of degree n — 1. 

System ([j]) has a unique singular point 0(0, 0) whose linear part of center 
type. 

Orbits of system ([[]) move around the origin with a constant angular velocity 
and the origin is a uniformly isochronous singular point. 

In Q the following problem was proposed: 

Problem 19.1. Identify systems w of odd degree which are O-symmetric 
(not necessarily quasi-homogeneous) having O as a (uniformly isochronous) cen¬ 
ter. 

We solve this problem for n = 5 and derive necessary and sufficient center 
conditions for the system 

x = y + x(ax 2 + bxy + cy 2 + dx 4 + ex 3 y + fx 2 y 2 + gxy 3 + hy 4 ), 
y = —x + y(ax 2 + bxy + cy 2 + dx 3 + ex 3 y + fx 2 y 2 + gxy 3 + hy 3 ), (||) 

a, b, c, d, e, f,g, h G R. 
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Theorem. The origin is a center of system & if and only if one of the 
following sets of conditions is satisfied: 


(i) a = b = c = 0,f = —3 (d + h)\ 

(ii) a = c = d= f = h = 0; 

(iii) a ^ 0, c = —a, f = 3 b(ae — bd)/(2a 2 ), 


(§ 


g = (2 a 2 bd + (2a 2 — b 2 ){bd — ae))/(2 a 3 ), 
h = (—2 a 2 d + b(bd — ae))/(2a 2 ). 


Proof. 


the comparison function j^] 


Necessity: 

To describe the behaviour of trajectories of @) near the origin we construct 


F{x, y ) = ( x 2 + y 2 )/2 + f 3 (x, y ) + / 4 (a:, y) + ..., 


where fk is a homogeneous polynomial of degree k whose derivative is 


dF 

dt 


D\(x A + y 4 ) + D 2 {x 6 + y 6 ) + D 3 (x 8 + y 8 ) + ... 


The number of the first coefficient Di other than zero defines the multiplicity 
of a complex focus and the sign of this coefficient defines stability of a focus; if 
Di = 0 for all i the origin is a center of (0. We refer to coefficients Di as the 
Poincare-Lyapunov constants. 

To find the Poincare-Lyapunov constants of a system x = p(x, y), y = q(x , y ) 
with a linear center we used computer algebra and wrote a Mathematica code 
that rests on the Poincare algorithm in |^|; see || for more details. 

PLconst[n_] := 


Module[{dF, ff, fF, x, y, pP, qQ, dD}, 
fF [2] := (x~2+y~2)/2; 

fF[i_] := Sum[ff[i-j, j]*x~(i-j)*y~j, {j, 0, i>]; 
pP[l] := y; 

pP[i_] := Sum[p[i-j, j](i-j)*y“j, {j, 0, i>]; 
qQ [1] := -x; 

qQ[i_] := Sum[q[i-j, j]*x~(i-j)*y~j, {j, 0, i>]; 
dF[k_] := (Sum[D [fF [i], x]*pP[k+l-i], {i, 2, k}] + 
Sum[D[fF[i], y]*qQ[k+l-i], {i, 2, k}])//Expand; 

Do [ 


Solve[Table[Coefficient[dF[k], x~(k-j) y~j], 

{j, 0, k>] 

==Table [0, {k+1}] , 

Table[ff[k-j, j], {j, 0, k>] 

]/.Rule->Set; 

Solve[Table[Coefficient[dF[k+1], x~(k+l-j)*y~j], 
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{j, 0, k+1}] 

==Flatten[{dD[k], Table[0, {k}], dD[k]>] 

&&ff [0, k+1]==0, 

Flatten[{Table[ff[k+l-j, j],{j, 0, k+1}], dD[k]}] 

]/.Rule->Set, 

{k, 3, 2n+l, 2}]; 

Table[Numerator[Together[dD[k]]],{k, 3, 2n+l, 2}] 

] 

The procedure PLconst [n] returns a list {Di,..., D n } of the Poincare- 
Lyapunov constants if we define the coefficients pij, qtj (2 < i + j < 2n + 1) in 
the Taylor series expansion of functions p(x, y ), q(x, y) beforehand. 

Using this procedure, we found the first four Poincare-Lyapunov constants 
of (|). 

Di = 2 (a + c), 

D 2 = —4ab — 46c + 3d + f + 3h, 

D 3 = 2(—85a 3 + 15a6 2 — 67a 2 c + 156 2 c + 61ac 2 + 43c 3 — 24 bd — 34ae— 

—22ce — 12b f — 50 ag — 38 eg — 48 bh), 

D a = 44600a 3 6 + 2736a6 3 + 84696a 2 6c + 2736 b 3 c + 47688 abc 2 + 75926c 3 - 
—37120a 2 d - 178 26 2 d - 32552acd - 2704c 2 d + 2364a6e + 12846ce- 
—2673de - 6120a 2 / - 2346 2 / - 3384ac/ + 792c 2 / - 891e/+ 
+6876abg + 50766c 5 - 3807 dg - 1269/ ff + 4720a 2 6 + 10986 2 6+ 
+31448ac6. + 19456c 2 6 - 2673 eh - 3807 gh. 

It is easy to verify that the equalities Di = 0,* = 1,2, 3,4 are equivalent to 
the following relations 

a T c = 0, 3d -} / T 3 h = 0, 

3ce — 6/ + 3 eg — 6 bh = 0, 2c 2 / — 3bcg + 3b 2 h = 0. 

If a = 0 then our simultaneous polynomial equations have two sets of solu¬ 
tions indicated in (i) and (ii). If a ^ 0 then, in view of the condition c = —a, 
we see that the other three equations constitute a non degenerate linear system 
for determining the variables f,g, h. The solution is given by (iii). 

The necessity part of the theorem is proved. 

Sufficiency: 

Case (i). System ([|) now takes the form 

x = y + x(dx A + ex 3 y + fx 2 y 2 + gxy 3 + hy A ) = y + xp±{x, y), ,ri. 

y = —x + y{dx A + ex 3 y + fx 2 y 2 + gxy 3 + hy 4 ) = x + yp 4 (x, y). 

This is a quasi-homogeneous system of degree 5 whose coefficients satisfy the 
equality / = —3(d + h ) which is the necessary and sufficient center condition in 
the case we study Q. 
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(I) 


Case (ii). System (^]) now takes the form 

x = y + x 2 y(b + ex 2 + gy 2 ), 
y = —x + xy 2 (b + ex 2 + gy 2 ). 

The planar differential system 

x=p(x,y),y = q(x,y) (f| 

is said to be reversible (in the sense of Zoliidek) if its orbits are symmetric with 
respect to a line passing through the origin. 

System d) is reversible if there is a linear transformation S : R 2 —> R 2 , 
sending a point ( x,y ) to the point (x',y') symmetric to ( x,y ) with respect 
to the line ax + /3y = 0 and satisfying the condition S(p(x,y),q(x,y)) = 
~{p(S(x,y)),q(S(x,y))). 

A more general condition of reversibility is as follows 

2 a/3(p(x, y)p(x',y') - q(x , y)q(x\ y'))+ 

(P 2 - a 2 )(p(x,y)q(x' ,y') + p(x',y')q(x,y)) = 0. 

It is well known that if system d) is reversible and has a linear center at the 
origin then the origin is a center of this system (see [[|] , for example). 

Obviously, system (|) is reversible because its trajectories are symmetric 
with respect to both coordinate axes. So, the origin is a center for system d). 

Case (iii). System (Q) now takes the form 

(2 a 3 )x = (2ci 3 )y + x(ax 2 + bxy — ay 2 ) 

(2 a 3 + 2 a 2 dx 2 — 2 abdxy + 2 a 2 exy + 2 a 2 dy 2 — b 2 dy 2 + abey 2 ), 

(2 a 3 )y = — (2a 3 )x + y{ax 2 + bxy — ay 2 ) 

(2a 3 + 2 a 2 dx 2 — 2abdxy + 2a 2 exy + 2 a 2 dy 2 — b 2 dy 2 + abey 2 ) 

H (B 

It turns out that system (Tj) is reversible. Its trajectories are symmetric 
with respect to each of the two perpendicular lines defined by the equation 
ax 2 + bxy — ay 2 = 0. The appropriate linear transformation S is given by each 
of the two matrices 

Su = ±( *,> + *)-■/.( ^ (»). 

This fact is confirmed by the straight calculations. We used Mathematica 
here. 

With the coordinate change x i— > xcosip + y sirup, y i— > — xsiny> + ycospi 
where the angle <p is defined from the condition atan 2 iy 2 + btaiup — a = 0, 
system (jrj) becomes as follows 

x = y + x 2 y(bi + exx 2 + giy 2 ), 
y = -x + xy 2 (bi + eix 2 + giy 2 ). 
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Hence the origin is a center for system (2j) in this case once again. 

The theorem is proved. 

2. It is known that isochronism of a center of a planar polynomial system is 
equivalent to existence of an analytic transversal system commuting with a given 
system in a neighbourhood of a center Js| ; observe that an arbitrary polynomial 
system with isochronous center not necessarily commutes with a polynomial 
system [| |. 

It is proved in || that if the systems 

x = p(x,y),y = q(x,y) ,g. 

x = r(x,y),y = s(x,y) ^ 

commute then y{x, y) = l/(p(x, y)s(x , y) — q(x , y)r(x, y)) is an integrating factor 
of both systems. 

Thereby if both commuting systems are polynomial then we can find the 
integrating Darboux factor for the given system and integrate the latter (about 
the method of Darboux and the relevant definitions see |J, f° r example). 

We now state the following fact which will be useful later. 

Considering (||), assume that 

p(x, y) = y + xR{x, y), q{x , y) = -x + yR(x, y ), 
r(x, y) = xQ{x , y), s{x , y) = yQ(x , y), 

where R(x,y),Q(x,y) are polynomials in a;,y. Then the algebraic curves x 2 + 
y 2 = 0, Q{x, y) = 0 are invariants for each of these systems. 

Indeed, it is immediately obvious that x 2 + y 2 = 0 is an invariant of both sys¬ 
tems with the cofactor 2 R(x, y) and 2 Q(x, y) respectively. The curve Q(x, y) = 0 
is an invariant of the second system with the cofactor xQ x (x, y) + yQ v (x, y). 

Because our systems commute the Lie bracket of the vector fields (p, q) and 
(r, s) vanishes and we have 

Px(x,y)r(x,y) +p y (x, y)s(x, y) - r x (x,y)p{x,y) - r v {x 1 y)q{x,y) = 0, 
or 

xQ{x, y)(R(x, y) + xR x (x , y)) + yQ(x, y)( 1 + xR y (x, y))- 
-p(x, y)(Q(x, y) + xQ x (x, y)) - xq(x, y)Q y (x, y) = 0, 
or 

x{Q x {x, y)p{x, y) + Q y (x, y)q(x, y)) = 

= (R(x, V ) + xR x (x, y))xQ(x , y)+ 

+(1 + xR y (x, y))yQ(x, y) - Q{x, y)p(x, y) = 

= {R{x, y ) + xR x (x, y))xQ(x , y)+ 

+(1 + xR y (x, y))yQ(x, y) - Q{x, y){y + xR(x, y)) = 

= x(xR x (x, y) + yR y (x , y))Q(x, y). 

We see that the curve Q(x, y) = 0 is an invariant with the cofactor xR x (x, y)+ 
yR v (x,y). 
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In this case p(x,y) = l/(Q(x, y)(x 2 +y 2 )) is an integrating Darboux factor. 


3. In each of the three cases we have found a non trivial polynomial system 
commuting with the respective system. 

In case (i) such a system is 

x = x(\ + ex 4 — 4dx 3 y + 4 hxy 3 — gy 4 ) = x(\ + q±(x, y)), 
y = y( 1 + ex 4 - Mx 3 y + 4 hxy 3 - gy 4 ) = y{ 1 + q 4 (x, y)). 


The function 

M(X,V) = (x 2 +y 2 )(l + q 4 (x,y))' 
is the integrating Darboux factor of (Q) and the function 


H{x,y) 


(x 2 + y 2 ) 2 
1 + q*(x,y) 


is the first rational integral of (||). 

The algebraic curves x 2 + y 2 = 0,1 + q 4 (x,y) = 0 are invariant curves for 

(I- 

According to || system (4) has a center of type B k , 1 < k < n — 1 whose 
boundary is a finite union of k unbounded open trajectories. Using (|^), in case 
(i) we can describe this boundary explicitly: 


1 

8 (o - ^(cos^.sin^)) 1 / 4 ’ 


where cq = maxr 0i2 7 r] (? 4 (cos ip, sin ip),x = g cos ip, y = gsin p. 

A straight analysis of this expression allows us to conclude that in our case 
a center may be of type B 2 or B 4 only. 

In case (ii) system (^!) commutes with the system 

x = (e — g)x + x{ex 2 + gy 2 ){b + ex 2 + gy 2 ), 

V = (e - g)y + y(ex 2 + gy 2 )(b + ex 2 + gy 2 ). 


This permits us to find an integrating Darboux factor 


y(x,y) = 


{x 2 + y 2 )(e — g + ( ex 2 + gy 2 ){b + ex 2 + gy 2 ))' 


The algebraic curves x 2 + y 2 = 0, e — g + (ex 2 + gy 2 )(b + ex 2 + gy 2 ) = 0 are 
invariant ones for system (|j). 

If b = 0 system (j|) is a system of the form ([!]) for which the condition 
/ = —3 (d + h) is obviously fulfilled. Then its first integral is 


H(x,y) 


(x 2 + y 2 ) 2 
1 + ex 4 — gy 4 
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If b ^ 0 then we may suppose that 6 = 1. The general case reduces to this 
by the change of variables x —> x/Vb, y —> y/Vb for 6 > 0 or x —> y/y/~^b, y —> 
xj\f—b, t —* —t for 6 < 0 . 

Then our system takes the form 


x = j/ + x 2 j/(l + ex 2 + 5 ?/ 2 ) = Xi(x,y), 
y = -x + xy 2 (l + ex 2 + gy 2 ) = Yi{x,y). 


( 0 ) 


The functions 

Ci = x 2 + y 2 , C 2 =e- g + ( ex 2 + gy 2 ) + (ex 2 + gy 2 ) 2 
are invariants for ( 0 ) with the cofactors 

Li = 2 xy{l + ex 2 + gy 2 ), L 2 = 2xy(l + 2 (ex 2 + gy 2 )). 
Moreover, if e ^ g the function 

2 . 2 

,y dt 


C 3 = exp( [ 
J o 


e — g + t + t 2 


is invariant with the cofactor L 3 = 2 xy. 

We have 2 Li — L 2 — L 3 = 0. Then the function 


H(x,y) = 


(x 2 + y 2 ) 2 


Cl _ 

C 2 C 3 (e — g + (ex 2 + gy 2 ) + (ex 2 + gy 2 ) 2 ) exp(/ f 


ex^+gy 2 dt 


is the first Darboux integral of (Hi- 
Let us remark that 


dt 


: arctan - 


1 + 2t 


J e - g + t + t 2 i/4(e- g) - 1 \/4(e - g) - 1 

for 4(e — 3 ) — 1 > 0 and 


/ 


e — 


dt 

g + t + t 2 


for 4(e — g) — 1 = 0 and 


2 

l + 2t 


/ 


e — 


dt 

g + t + t 2 


1 l + 2t- x /l-4(e-g) 

-\/l — 4(e — g) 1 + 2 1 + y/l — 4(e — g) 


for 4(e — g) — 1 < 0. 

If e = g the function 


1+* 2 

C 3 = exp x ' 2 +y ' 2 
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is invariant with the cofactor for © £3 = -2exy. 
We have 2Li — L 2 + ^£3 = 0. Then the function 


H(x,y) 


cl 

C 2 C 3 1/( 


(x 2 + y 2 ) 

(1 + ex 2 + ey 2 ) 


1 + a; 2 , 

6Xp( 2 1 9 ) 

e ar + 


is the first Darboux integral of for e = g. 

Since (0) has a unique finite singular point at the origin, the phase portraits 
are obtained by studying the points at infinity. A standard inspection of the 
location and types of such points on the equator of the Poincare sphere allows 
us to conclude that (0 has phase portraits of two types only: a center is of 
type B 2 when eg > 0 or of type B 4 when eg < 0. 

In case (iii) a commuting system and first integral may be found on consid¬ 
ering that system (Q) is equivalent to (0. 

Observe that for d = e = 0 system (Id) is a quasi-homogeneous O-symmetric 
cubic system of the form 


x = y + x(a x 2 + bxy — ay 2 ), 
y = —x + y(ax 2 + bxy — ay 2 ). 


It commutes with the system 


x = x + x(bx 2 — 2 axy), 
V = y + y(bx 2 - 2 axy), 


and has the first integral 


H{x,y) 


x 2 + y 2 

1 + bx 2 — 2 axy 


Summarizing we conclude that the system under consideration has phase 
portraits of two types only: a center is of type B 2 or of type B 4 . 
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